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Abstract 

We show that the evolution equations in QCD predict geometric scahng for quark and 
gluon distribution functions in a large kinematical window, which extends above the satu- 
ration scale up to momenta of order 100 GeV^. For < Ql, with Qs the saturation 
momentum, this is the scaling predicted by the Colour Glass Condensate and by phenomeno- 
logical saturation models. For 1 ^ ln{Q^/Ql) <C \n{Ql / Aqq^)) , we show that the solution to 
the BFKL equation shows approximate scaling, with the scale set by Qs- At larger Q^, this 
solution does not scale any longer. We argue that for the intermediate values of where 
we find scaling, the BFKL rather than the double logarithmic approximation to the DGLAP 
equation properly describes the dynamics. We consider both fixed and running couplings, 
with the scale for running set by the saturation momentum. The anomalous dimension which 
characterizes the approach of the gluon distribution function towards saturation is found to 
be close to, but lower than, one half. 



1 



1 Introduction 



At very high energies, the contribution to the hadronic wavefunction which dominates hadronic 
processes corresponds to a state of very high partonic (mostly gluonic) density. This high density 
matter is beUeved to reach a saturation regime, and become a Color Glass Condensate ||l|, ^. 
(See also for a recent review and more references.) Recent phenomenological analysis of the 
data from HERA and the RHIC show results qualitatively and semi-quantitatively consistent 
with this picture Q. The results, while not entirely compelling, provide a strong motivation for 
the further theoretical investigation of this novel form of high density matter. 

In a very important paper Stasto, Golec-Biernat and Kwiecihski have shown that the 
HERA data on deep inelastic scattering (DIS) at low x ||6|, which are a priori functions of two 
independent variables — the photon virtuality and the Bjorken variable consistent 
with scaling in terms of the variable 

T = Q^Rl{x) (1.1) 

where Rl{x) = {x/xq)^/QI with A = 0.3 - 0.4, Qq = I GeV, and xq ~ 3 x 10"'' in order to 
fit the data. In particular, the data for the virtual photon total cross section at x < 0.01 and 
< 400 GeV^ are consistent with being only a function of T. Note that, because a^*p = 
47r^aEM-^^2(ic, (5^)/Q^, this implies that F2{x,Q'^)/Q'^ is a function only of T in the indicated 
kinematical range. 

At first sight, it is tempting to interpret this scaling as a consequence of the existence of the 
Color Glass Condensate: Indeed, it has been argued that in the saturated region the various 
parton densities are functions only of Q"^ /QHx), with Qg the saturation momentum, which is 
proportional to the gluon density, and thus increases rapidly with 1/x. Although Qg has not 
been fully determined from first principles, its estimates based on the quantum evolution of the 
Color Glass Condensate suggest indeed a power-like increase: ~ QqX^''^ [0, ^, ^, |l^, pT]] . 
(This will be further discussed in this paper.) The fundamental problem with this argument is 
that it is valid only for less than or of the order of the saturation momentum, which is at 
most several GeV^, while the fit of Ref. |5| extends up to of the order of several hundred 
GeV2. 

The main purpose of this paper is to show that the scaling region for the various distribution 
functions is in fact much larger than the saturation region. We shall find geometric scaling for 
all momenta such that the following inequality is satisfied: 

ln(QVQ^)«ln(Q?/A^CD)- (1-2) 

For ~ 2 GeV and Aqcd ~ 200 MeV, the upper scale on Qj/A^^^.^ ~ 400 GeV^. While 
at soft momenta ^ this scaling is an expected consequence of saturation, at high momenta 
1 < ln{Q^/Qf) <C ln(Qg/AQQj-,) it rather corresponds to a regime where parton densities are 
small, and linear evolution equations apply. And indeed we shall see that the extended scaling 
at Q2 > q2 

arises from solutions to the BFKL equation jl^, which is the appropriate limit of 
the general non-linear evolution equations [2,13-18] in the kinematical range of interest. 

Note that, in this paper, we shall use the BFKL equation as an effective equation, valid in 
some range of whose increase in is correlated with the decrease in Bjorken x (because of the 
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constraint > Q'l{x) ~ QgX^'^). Thus, conceptual difficulties of BFKL, like the exponential 
increase of its solution at high energy, and its "infrared diffusion", will be of no concern for us. 
In fact, an essential step in our analysis will be to clarify the kinematical ranges in which the 
solution to the BFKL equation (in the saddle point approximation) is of the "genuine BFKL 
type" (i.e., it is close to the saddle point which governs the limit x — > with fixed), or of 
the "double logarithmic type" (i.e., close to the saddle point which describes the limit — > oo 
with X fixed). As well known, in the latter limit the BFKL solution becomes equivalent to the 



double logarithmic approximation (DLA) to the DGLAP equation |19]. But this limit turns out 
not to be relevant for the geometric scaling: Indeed, the window for extended geometric scaling 
that we shall find is actually in the range controlled by the BFKL saddle point. 

It is quite remarkable that even at high momenta > Q^, it is still the saturation momen- 
tum Qsix) which sets the scale in the scaling variable (i.e., which plays the role of l/Ro{x) in 



eq. (LI)). This is so because the solution to the BFKL equation must properly match when 
Q1 on the corresponding solution at saturation. 

In brief, our strategy in this paper will be as follows: We start with a solution to the BFKL 
equation, which in general (i.e., for arbitrary Q"^) does not appear to scale. By extrapolating this 
solution down to momenta ~ Q'^g{x)^ we shall estimate the saturation scale Qs{x) from the 
matching condition with the solution at saturation. The latter is not precisely known around 
Qsi but for the present purposes such a precise solution is not really needed: All that we need is 
a saturation criterion which gives us the value of the quantity of interest at = Ql{x). This 
criterion is simply that the gluon distribution function itself becomes of order l/ag. Such a 
criterion follows from the general analysis of the non-linear effects due to high parton densities 
[1-3,7-11,13-18,20,21]. By expanding the BFKL solution around = Ql{x), we shall find 
geometric scaling in the momentum range where the first term in this expansion is a good 
approximation, namely for 1 < \n{Q'^ /Ql) <C \d.{Q1 / K^qy)) ■ This conclusion holds both for 
a fixed coupling and for a running coupling (with the scale for running set by the saturation 
momentum), although the saturation scales turn out to be different in the two cases. 

Note that, although the saturation is crucial for the arguments in this paper, the emergence 
of an extended scaling window above the saturation scale is not an automatic consequence of the 
saturation. Rather, this requires some non-trivial properties of the linear evolution equations as 
well. To see this, consider some dimensionless functions f{x,Q'^) (e.g., the gluon distribution) 
which at high obeys a linear evolution equation (BFKL or DGLAP), but at ~ Ql{x) 
satisfies a saturation condition of the type (after an appropriate normalization): 

f{x,Q^ = Ql{x)) = 1. (1.3) 

As a consequence of this condition, and without loss of generality, the function / near = Ql{x) 
(i.e., for I ln((5^/Qs)| <C ln{Q'^ / Aqq^))) , can be approximated as: 

where Xs{x) = {din f /dlnQ'^)\Q2(^^y Assume that this limited expansion can be extended up to 
values which are high enough for the non-linear effects to be subleading. Then, the exponent 
\s{x) can be determined by solving the BFKL equation with the saturation boundary condition 
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( |1.3| ). A priori, the result Xs{x) can be a non-trivial function of x. If this happens, then eq. (|TT 
shows no geometric scaling. However, in our subsequent analysis, we shall discover that for 
solutions to the BFKL equation. As is a number independent of x — as a consequence of the 
scale-invariance of the BFKL kernel — , so that /(x, Q^) shows scaling in the window of validity 
of the expansion ( |1.4| ). The difference between the actual value of the exponent and its 
"naive" value expected from perturbative considerations at large will be referred to as the 
"anomalous dimension" . 

The saturation criterion that we shall actually use refers to the scattering of a small dipole 
off a hadronic target: 1/Qs is the critical transverse size of the dipole at which the scattering 



amplitude becomes of order one ("black disk limit") |2^, 16, 23, ^. This criterion, which is 
equivalent (at least, for our present purposes) with the more standard condition on the gluon 
distribution |l], |, |, |l|, |l|, |2|, is simpler to use in applications to deep melastic scattering. 



since more directly related to the corresponding cross-section at small x ||2^, 0, |2^, 

cTT,L{r,Q^) = J^'^^J ^^V|*T,L(^,r±;Q2)|2^(r,rx). (1.5) 

In this equation, (Tt,l{'^iQ'^) is the cross-section for the scattering of a virtual transverse (T) 
or longitudinal (L) photon off the hadron, at relative rapidity r = ln(l/x) ~ Ins. (Note that, 
from now on, we use r to indicate the dependence upon the total center-of-mass energy squared 
s.) Furthermore, is (minus) the photon virtuality, and ^'t,l(-2, Q^) is the light-cone 
wavefunction for the photon splitting into a qq pair (the "dipole") with transverse size r_|_ and 
a fraction z of the photon's longitudinal momentum carried by the quark. Finally, (t(t, r_L) is 
the dipole-hadron cross-section, which can be computed in the eikonal approximation as: 

a(T,rx) = 2J d^b^{l-Sr{x±,y±)), Sr{x±,y±) = ^ {ti{vHx^)V{y±)))r, (1.6) 

with r_L = x_L — y± (the quark is at x±, and the antiquark at y±), and the impact parameter 
b± = {x± + y±)/2. In eq. ( |1.6| ), V"^ (or V) is a Wilson line along the straight line trajectory of 
the quark (or the antiquark), that is, a path ordered exponential of the color field created in 
the hadron by color source at rapidities r' < r. The brackets in the definition of the 5'-matrix 
element 5,- {x± , y± ) refer to the average over all configurations of these colour sources Q . 

A small dipole is weakly interacting: Sr{r±) ~ 1 for r± <^ 1/Qs{t). A large dipole, on the 
other hand, is strongly absorbed: Sr{r±) <C 1 for r± ^ 1/Qs{t). This is so because of the large 
density of the saturated gluons in the hadron wavefunction. At momenta ^ Q'^{t), there is 
only one intrinsic scale in the problem, the saturation momentum Qs{t) itself (we mean this for 
a fixed impact parameter and a fixed coupling). So, all physical quantities should be expressed 
as a dimensionless function of Q"^ /Ql times some power of giving the right dimension. These 
general properties expected for a{T,r±_) have been incorporated in a simple "saturation" model 
by Golec-Biernat and Wiisthoff ||2^ {1/Rq{x) plays the role of the saturation scale) : 



.(.,.x) = .o(l-e-i/".^«).,„.(j^), (LT) 

which appears to give a reasonable description of the HERA data at small x. By inserting 
this expression in eq. ( |l.5D , one can check that the scaling property is transmitted to the DIS 
cross-section: c^*p(t, Q^) = a^*p{T) with T of eq. ( |1.1| ) ||5|. 
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shows exact scaling for all distances r±, and not only in the 
In reality, however, we know this scaling to be violated at 



Written as it stands, eq. ([T 
saturation regime r± S> Ro{x) 
sufficiently small distances r± <^ Roi^), where eq. (1/7) should be replaced by 

2 vr^a^ xG{x, 1/rj 



4 26, 



o"(T,r_L) = ao 1 - exp 



N, 




where xG{x,l/r'j_) is the gluon distribution function evaluated at = » Q^. At such 

high momenta, xG{x, Q^) is a solution to some linear evolution equation, usually taken to be the 
DGLAP equation, which, at least for sufficiently large Q^, has no scaling behaviour (see Sect. 
3 below). It will be our main objective in what follows to establish up to which momenta 
the scaling property holds for the dipole cross-section, and thus for the quark distribution in the 
hadron. To characterize this property also for the gluon distribution, we shall find it convenient 
to introduce a definition of the latter in terms of the dipole scattering amplitude (see eq. (|0 
below) . 



and eqs. (|2.7D -(p 

To simplify the problem, we shall assume that the dependence upon the impact parameter 
is not essential for the present purposes, so we can treat the hadron as being homogeneous in 
the transverse planeQ. This implies 5*1- (x_l , y_L ) = Sr{r_i_), with r_|_ the size of the dipole, and the 
scaling properties of cj^*p(r, Q"^) are directly related to corresponding properties of Sr{r_[_), that 
we shall study. 

The function ST-{r±_) will be obtained by solving an appropriate evolution equation in r. As 
already mentioned, we shall be mainly concerned with solutions to the BFKL equation, but this 
equation will be viewed as the linear limit of more general, non-linear, evolution equations. This 
is important, as it will allow us to keep trace of the non-linear effects via the boundary condition 
at ~ Q1{t) (the "saturation condition"). 

It is possible to write down a formal evolution equation for Sr{x^,y^) without specifying 
the average on the target hadron in the definition (|1.6|) JTsI, 16, |l8i. In general, such an equation 



is only the starting point of a hierarchy of coupled evolution equations for the correlators of 
Wilson line operators |15], which can be encoded in a single, functional evolution equation [ p^ . 
But a closed equation for >S't-(x_l, yx) can still be obtained in the large N^. limit. This has been 
derived by Kovchegov [l^ within the Mueller's dipole model [28|, and follows also from the 
general evolution equations by Balitsky |15| by taking the large Nc limit. (See also ||l^ for a 



different derivation.) We shall refer to this as the Balitsky-Kovchegov (BK) equation. 

Alternatively, an effective theory can be constructed for the small-x component of the hadron 
wavefunction [||, ^, ^,|2^, 29, 3^, thus allowing a direct and explicit calculation of the average in 
eq. (^]^) . In this effective theory, the high density gluon configurations at small x are treated as a 
"Colour Glass", i.e., as the classical colour field generated by colour sources at rapidities r' < r = 
ln(l/x), which are "frozen" in some random configuration. In this picture the evolution is viewed 
as a renormalization group operation in which layers of quantum fluctuations are successively 
integrated out to generate the sources |30, ^. This leads to a functional renormalization group 
equation (RGE) for the probability distribution of the colour sources. Remarkably, this RGE is 
equivalent to the Wilson line approach of Refs. |]l^, |l^, |l^: it generates the same evolution 

^This assumption is implicit in eqs. (1.7) and (1.8), where ctq = 2 J d?b± = 2ttR^, with R the hadron radius. 
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equations for the correlation functions of Wilson lines. At high transverse momenta Q^, the non- 
linear effects are weak and can be expanded out. To lowest order in this expansion, the general 
RGE (and also the BK equation) reduces to the BFKL equation jl^, as already mentioned. 
The outline of this paper is as follows: 

In the second section, we rely on the BK equation to explore the scaling properties of 
distribution functions in the deeply saturated region where the gluon density is of order 1 joLg . By 
assuming that scaling solutions exist, we determine the dependence of the saturation momentum 
upon T for both the fixed coupling and the running coupling cases. Throughout this paper, the 
scale for the running of is always the saturation momentum. 

In the third section, we consider the intermediate region where the evolution equation 
linearizes and reduces to the BFKL equation. We construct approximate solutions for this 
equation in the saddle point approximation and show that the solution which is relevant for 
the approach to saturation and the extended scaling is closer to the standard BFKL solution, 
rather than to the double-logarithmic DGLAP-like solution. We determine the saturation mo- 
mentum, compute the anomalous dimension of the distribution function, and show the solution 
has geometric scaling for both the cases of running and fixed coupling. 

The last section contains a summary of our results, and some discussion. 



2 Scaling properties of the Balitsky-Kovchegov equation 



In this section, we shall briefly discuss the BK equation |15| , 16], which is a non-linear equation 
for the evolution of the dipole-hadron scattering amplitude with r = ln(l/x). This equation will 
be used here just for qualitative arguments, which by themselves are valid for any although 
the BK equation holds, strictly speaking, only in the large limit. (The same arguments at 
finite Nc would require the full formalism in Refs. []l5| , 0, |2|, |3|.) Specifically, the following 
properties of the BK equation will be useful for what follows: First, it covers both the linear 
and the non-linear regimes of the quantum evolution (according to the value of second, 
it reduces to the BFKL (and eventually DLA) equation in the linear regime at high Q^, and, 
third, it demonstrates the relation between geometric scaling and saturation at low . 

In fact, we shall find that the structure of the BK equation is consistent with scaling solutions, 
although the kind of global arguments that we shall use cannot decide if such solutions actually 
exist or not. This will open the discussion of scaling solutions outside the saturation regime, to 
be pursued in the framework of linear evolution equations in Sect. 3. The scaling properties of 
the Balitsky-Kovchegov equation have been also investigated numerically, in Refs. |^, |ll|, 

The BK equation is most succinctly written as an equation for the S'-matrix element ^^(rx) 



for dipole-hadron scattering, eq. (|l.6|) (below, = Ncag/i^) 



dr J 271 (x_L - 2:_L)2(y_L - 

X {Sr{x^ - y±) - Sr{x± - z±)Sr{z± - y±fj . (2.1) 

Given Sr{r±), the dipole-hadron scattering amplitude A/'r(r_L) is obtained as: 

K{r±) = 1- Srir^). (2.2) 
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For a small dipole, r± <C 1/Qs{t), Nr{ri_) is small as well, A/'r(?"±) ^ 1, and one can linearize 
eq. (|2.1| ) with respect to Nr{T±): one then obtains the BFKL equation. This can be easily seen 
in the momentum space form of the BK equation. (The BFKL equation in coordinate space will 
be discussed in the next section.) Specifically, if one defines (/3t(^±) by 



rl J (27r) 
one finds 

The small-dipole condition A/'r(?^±) <C 1 in coordinate space corresponds to the condition 
^t{'1±)/<1± <C 1 in momentum space, which is satisfied provided q'j_ is large enough, q'j_ ^> Q1{t). 
In this regime, the last term, quadratic in ^t{q±)/q'±, can be neglected in the r.h.s. of eq. ( |2.4| ), 
which then reduces to the BFKL equation For momenta q'j_ which are even larger, this 

becomes equivalent to the double-log limit of the DGLAP equation (see Sect. 3 below). 

On the other hand, non-linear effects become crucial at low momenta q'j_ ^ Q'^{t), that is, for 
a large dipole which is strongly absorbed by the hadron: Sr{r±) <^ 1 for r± ^ 1/Qs{t). In this 
regime, the right hand side of eq. (2.1) can be simplified by neglecting the term quadratic in Sr 



this is appropriate since the dominant contribution (in the sense of the leading log) comes from 
z± satisfying 1/Qs(r) <C |2± — <C r_L (or a similar condition on \z± — y±\). The equation 
then reduces to 

-^Srir^) ^ -aAn[rlQl{T)) 5.(r^), (2.5) 

for which a solution will be written down shortly. (This requires the r dependence of the 
saturation scale Q'^{t).) 

Still in the non-linear regime, it is interesting to determine the limiting form of the function 
(/7t-(9±) at low q±. By solving eq. (|2.4| ) in this regime, or, more directly, by noticing that 
■N'rifi.) — 1 for r_L > \/Qs{t) in eq. ( ^.31) , one deduces that: 

ipr{q± < Qs) ^ q± / = 7rq±ln — —. (2.6) 

It can be checked that this is indeed a solution to eq. ( ^.41) at low q'^ and to leading-log accuracy 



At this stage, one may notice that the properties of the function iPt{q±) introduced in 
eq. ( |2.3D are very similar to those expected for the "unintegrated gluon distribution", i.e., the 
density of gluons in the transverse phase-space. Indeed, at large momenta ^PriQ±) satisfies the 
BFKL equation and is truly proportional to the gluon distribution, as manifest on eq. ( |1.8| ). At 
low momenta, q'j_ <C Q^j it has the same behaviour in q±, eq. (^), as the gluon density in the 
saturation regime [0, ^] . It is thus quite reasonable to define the gluon phase space distribution 
function via the Fourier transform of the dipole scattering amplitude (with Cp = {N^ — l)/2Nc): 

(9 ^2 ^ ^r{k±,b±) 
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where we have reintroduced the impact parameter dependence, for more generahty (so, 
^T{k_\_-,h_i) /k\_ is actuaUy the Wigner transform oiJ\fr{r_\_, h_[_)/r\)^ and the overall normalization 



follows by comparison with eq. (1.8). This gives a gluon distribution: 



dr TT^Us J {2-kY J J rj_ 

The definition ( |2.7| ) should be compared to the canonical definition of the gluon density, which 
involves the expectation value of the gluon occupation number in the infinite momentum frame 
and in the light-cone gauge |, |,|. In the hnear regime at high Q^, these two definitions are 
equivalent, and provide both the standard gluon distribution which measures (Bjorken) scaling 
violation in F2. But in the non- linear regime at ^ Q^, there is no simple relationship between 



F2 and the canonical gluon density. By contrast, the definition (2/7) has the advantage that it 
appears linearly in the expression for F2, so it is related to a directly measurable quantity even 
at saturation. In particular, all the scaling properties that we shall establish later for the dipole 
scattering amplitude MT{r±) will immediately translate to the gluon distribution defined as in 
eqs. ( p.7[ )-( [2^ ), and we expect them to hold for the canonical gluon distribution as well. 

Note the factor of l/a^ in the relation (|2.8| ) between Nt{t±) and xG{x, Q^). Via the running 
of as, this will be a source of geometric scaling violation in xG{x,Q'^), as we shall show later 



that A/'r(r_L) shows scaling in the range specified by eq. (1.2). 

Eqs. ( p.5| ) and ( |2.6D confirm already our expectation that, in the saturation regime, the 
dipole scattering amplitude should exhibit geometric scaling: the solution Sr{r±) to eq. ( p.5| ) is 
clearly a function of r'j_Q'^{T) alone, and, similarly, the dimensionless quantity ^t{q±)/q± is only 
a function of q^/Qlir) when q'j_ <C (5s(t). Encouraged by this observation, let us investigate 
the scaling properties of the solution to the BK equation in more generality: 

1) Fixed coupling case: 

We consider first the case where the coupling as in the r.h.s. of eq. ( |2.1[ ) is fixed, and search 
for a solution Sr{r±) to this equation in the scaling form: 

Srir^)^l-m, e^ln^;^^. (2.9) 

Assuming such a scaling, the r dependence of the saturation scale is then fixed by the equation. 
To see this, note first that for any function /(^), 



fio = - (g^i^Q'c^)) /'(o, ri^m = -/'(e), 



dr 

with /'{£,) = df{^)/d^. Then, integrating the BK equation (|2.1|) over r_\_ = x_\_ — y± (after first 
dividing it by r'j_), one gets 

/ dV^^5r(r±) =vr(5.(cx,) -5.(0))^lnQ2(^) = _^l_inQlir), 

where we have used the boundary conditions S'r(O) = 1 and S't-(c>o) = 0. Therefore, 

^lnQl{T) = cas, (2.10) 
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where c is given by 

d?r^(Pz_[_ 



27r2 (r_|_ — ] 



Srir±) - Sr{z±)Sr{r± - ZJ_ 



(2.11) 



If Sr{r±) is a scaling solution, then the r.h.s. of eq. ( 2.11| ) is a constant independent of r. 
This follows from the scale invariance of the integrand: The function (2^) depends upon r 
only via the scale Qs{t) within the scaling variable; thus, by changing variables according to 
u\ = r\Qs{T) and v\ = z]_Qs{t), all the r dependence goes away. More explicitly, the integral 



in eq. (2.11) can be written only in terms of the scaling variable ^ = In l/rj_Qj{T): 



HO - m 

le-«' -e-?l 



+ 



-2? 



(2.12) 



Therefore, the r.h.s. of eq. ( ^.10 ) is independent of r, which implies that the saturation scale 
grows exponentially with r : 



A2 e^°^^, 



(2.13) 



with A fixed by the initial condition (typically, A ~ Aqcd)- With this saturation scale, eq. ( |2.5| ) 
can be easily integrated, with the following result which shows how the black disk limit (Sr = 0) 
is approached when S> l/Qs |3| : 



Sr{r±) oc exp |-|-- 



(2.14) 



In order to search for scaling solutions, it is preferable to rewrite the original BK equation as 
an equation for the scaling function <&(^). It turns out that the equation takes a simpler form 
when written in momentum space. That is, the scaling Ansatz (|2.9| ) is reformulated at the level 
of the BK equation in momentum space ( |2.4[) . In view of eq. ( ^.31) , it is natural to write: 



*(C) 



C = In- 



1l 



The function ^'(C) is related to ^{C) via the following relation, which follows from eq. (|2 



(2.15) 



-oo 47r 



3^-«Jo(e(f-«)/2)M/(C) 



where Jq is a Bessel function. By inserting the Ansatz ( 2.15| ) in eq. (| 
eq. ( 2.10 ), one obtains the following equation for the scaling function ^'(C): 



c|^(C) 



dC' 



■ef'^'(C') -e^^'(C) 



=C' -eCl 



+ 



e^vI/(C) 
V4e2C' + e2C J 2tt 



M/2(C) 



(2.16) 
and using also 

(2.17) 



This equation, together with the definition ( 2.12| ) of c, and the relation ( ^.16| ) between ^(C) and 
<I>(^), form a system of a coupled equations which in principle determine the scaling solution 
and the associated coefficient c. 



2) Running coupling case: 
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One should stress that the running of the QCD couphng is a higher order effect, which 
so far has not been included in the derivation of the non-linear evolution equations from first 
principles. Thus, treating in the BK equation ( |2.1D as a running coupling "constant" is just 
a phenomenological way to incorporate some (potentially large) higher order corrections, and 
suffers from ambiguities. Here, we shall assume the one-loop-like running 



bo 



UNr 



with the scale chosen as the saturation momentum: 



(2.18) 

(r) . This is physically acceptable 
since Qs is the typical momentum of the gluons in the saturated regime. Moreover, numerical 
solutions to the fixed coupling BK equation show that the dipole scattering amplitude is peaked 
around the saturation scale [11|. (Other possibilities for running, e.g., = l/r'j_, with r± the 
size of the dipole, will be left for a later study |]3^ .) 

Since the running coupling ( 2.1§| ) involves also the QCD scale Aqcdj this is not a single scale 
problem any longer. It is therefore a nontrivial question in general whether a scaling solution 
exists with running coupling. But at least for very low q'j_, the solution (Pt{q±) in eq. ( |2.6D is 
manifestly a scaling solution even for a running coupling, since independent of ag ! (Together 
with eq. (2.7), this reflects the fact that the gluon density at saturation is of order l/a^.) 
Thus, even in this case, it is legitimate to look for a scaling solution, of the form dz!" 



By 



the same arguments as above, this leads us to the differential equation 

casiQlir)), 



(2.19) 



where c is again constant, and takes the same value as in the fixed coupling case, since determined 
(together with the scaling function '&(i^)) by the same coupled equations ([2T2D , (|1|) and (|T^). 
As compared to the fixed coupling case, the growth of the saturation scale becomes somewhat 
milder (tq is an arbitrary constant) 



Qs(r) = A^CD eV2^oc(r+.o). (2.20) 

It would be interesting to clarify if the system of equations (2.12), (2.16) and ( 2. IT] ) has 
indeed solutions, i.e., if exact scaling solutions to the BK equation really exist. But as we 
shall argue in the next section, the physically relevant solutions necessarily violate scaling at 
sufficiently small r±. This is so since, when r±_ <^ 1/Qs(t), the BK equation linearizes, so its 
solution should match on the corresponding solution to the BFKL (or DLA) equation, which 
shows scaling only in a limited range of r± below the saturation length 1/Qs(r). 



3 The BFKL equation in the context of saturation 

We have seen in the previous section that, when the size of the dipole is small compared to 
the saturation scale, r± <C 1/Qs{t), the general evolution equation can be linearized in the 
scattering amplitude A/V {r± ) , and then it reduces to the (coordinate space form of the) BFKL 
equation [||, H : 
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The solutions to this equation have been extensively studied (see, e.g., |33] for an approach 
similar to ours), but they will be reconsidered here in the context of saturation, which requires 
the function A/'r(?'±) to satisfy the boundary condition A/'t-(?'±) ~ 1 for r± ~ 1/Qs(r). This is 
automatically satisfied by the solution to the non-linear BK equation, but in the framework of 
the linear BFKL equation it becomes a non-trivial condition which determines the saturation 
scale. (A similar strategy to determine Qs in the context of the BFKL equation has been 
previously used by A. Mueller [Q].) 

To understand this boundary condition, recall that, at saturation, Mr{r±) = 1 — Sr{r±) is 



a scaling function (see, e.g., eq. ( 2.14]) ): A/'r(?'±) = fifj_Qi{T)), and thus becomes a constant k 



when r±Qs{T) = 1. This constant is a number of order one (although strictly smaller: k < 1), 
since ST-{r±) <C 1 at saturation. The precise value of k is a matter of convention — it defines 
what we mean exactly by "the saturation scale" — , but this will not matter for what follows. 

3.1 The solution to the BFKL equation revisited 

In this subsection we shall construct approximate solutions to the BFKL equation ( |3.lD for the 
case of a fixed coupling constant a^. Although some of the results are quite standard (see, e.g., 
|l33|| ), we prefer to go through their derivation in some detail, in order to clarify the range of 
validity of the various approximations. This will be important for the discussion of extended 
scaling in the next subsection. Moreover, the techniques that we introduce here will be also 
useful later. 

The approximations that we shall perform rely on the following inequalities: 

n2 

OsT > 1, and r = ln-f75- > 1, (3.2) 

where A is some reference scale of order Aqcd ; and = l/rj_. The conditions (^]^) express the 
fact that we consider a perturbative regime at small x and large Q^. In addition, the quantities r 
and oisT are not free to vary independently; rather, they are constrained by ^ Q^{t) (which 
ensures that we are in a linear regime), which requires (cf. eq. ( |2.13| )): 

r = In --5- > In — -5— = ca^r , (3.3) 

with the coefficient c to be determined in Sect. 3.2. 

Note first that eq. ( p.lj ) has the same structure in coordinate space as the usual BFKL 
equation in momentum space (i.e., the linear part of eq. (|2.4| )), so it can be solved via similar 
techniques. For the present purposes, it is convenient to use the Mellin transform with respect 
to the transverse coordinates: 

where £^ = 1/A^ and r'j_ <C so that the contour is taken on the left of all the singularities 
of the integrand in the half plane Re A > 0. Since the BFKL equation is invariant under scale 
transformations, the ensuing equation for Xr(A) is local in A : 

— Xr(A) = a.{2V^(l) - V(A) - V(l - A)}xr(A), (3.5) 
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is the di-gamma function), and has the obvious solution 

X^{X) = ea=-{2V'(i)-^(A)-V'(i-A)} ^^^x). (3.6) 

The initial condition Xo('^) is not important for what follows (see below), so it is left unspecified. 
In order to return to coordinate space, we need to perform the integral 

[ ^ ^aM2^{i)-^{X)-^{i-x)}^^^x) = I ^e^(^'^'-) (3.7) 
Jc 27ri Jc 2-711 

where r = ln(r^/^'^) is negative and large in the range of interest. (That is, in coordinate space, 
the second condition ( p.2| ) is rewritten as (— r) ^1.) The function ( |3.6| ) has essential singularities 
at all the positive integers A > 1, so we can choose the contour as C = {A = a+iv, — oo < u < oo}, 
with < a < 1. 

The integral (|3.7| ) will be evaluated in the saddle point approximation, which is a good 
approximation when r and a^r are both large. (The corrections to it are suppressed by powers 
of 1/r or l/oisT.) To the same accuracy, we can ignore the initial condition Xo('^)) since its 
contribution to the function -F(A, r, r) in the exponent is not enhanced by either r or a^r. (We 
implicitly assume here that Xo(-^) is not rapidly varying in the range of A of interest.) 

We thus obtain: 

i-oo 27r ^2ttF"{Xo) ' ^ ' ' 

where, from now on, 

F{\) = rA + asT {2^(1) - V(A) - ^-(1 - A)} = F(A, r, r), (3.9) 
and Ao satisfies the saddle point equation: 
(9F(A,r, r) 



= 0, Ao = Ao(r,r). (3.10) 

Ao 

In eq. ( |3.S|) , we have also included the contribution of the Gaussian fluctuations around the 
saddle point. That is, we have taken the contour C = {A = Aq + iv, — oo < v < oo}, and we 
have expanded in powers of v to quadratic order before integrating over v. 

To visualize the solution to eq. ( 3.10| ), it is useful to keep in mind that 2^/^(1) — -(/^(A) — '0(1 — A) 



is a convex function of A with its minimum at A = 1/2 and simple poles at A = and A = 1. A 
good approximation to this function in the range < A < 1 is given by 

2-0(1) - V(A) - V(l - A) « i + + 41n 2 - 4. (3.11) 

A 1 — A 

Thus, there is an unique saddle point Ao in the region < A < 1, whose position moves between 
and 1 depending upon the value of the ratio r/a^r. There are three limiting cases of interest: 

(A) When rjoLsT is positive and large enough: the saddle point is close to A = 0. 

(B) When rjagT is small ~ 0: the saddle point is close to A = 1/2. 

(C) When r/agT is strongly negative: the saddle point is close to A = 1. 
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More precise boundaries between these various cases wiU be specified below. 

The first case is relevant when we consider the BFKL equation in momentum space, where 
r = ln(A;^/A^) is always positive in the range of interest. Similarly, case (C) applies only to the 
coordinate space, where r = ln(r^/^^) is negative. These two cases, (A) and (C), correspond 
to the double logarithmic approximation (DLA), which describes the leading behaviour of the 
solution at large k\ (or small r^) for fixed, but large, r; this limit is common to the BFKL and 
DGLAP equations. Case (B), on the other hand, applies to both the momentum space (r > 0) 
and the coordinate space (r < 0). In standard analyses of the BFKL equation, this is the case 
which describes the high energy limit (r ^ oo at fixed r) ||3^. Here, we shall never be truly in 
this limit, because of the condition ( |3.3D which, strictly speaking, implies that rjoLgT is never 
small ! If case (B) is nevertheless relevant for us here, it is because the true saddle point remains 
close to A = 1/2 up to relatively large values of r/SsT, which, as we shall see, leaves enough 
space for the condition (^.31) to be satisfied. 



(A) Assuming that Aq ^ 1, we find the saddle point 



Ao - \/^, (3.12) 



which after insertion in eq. (3.8) and performing the Gaussian integral there, leads to 



I ~ e^v'^J^O— . (3.13) 

When r = ln(A;^/A^), this gives the well-known "DLA solution" in the momentum space, which 
coincides with the solution to the DGLAP equation in the double-log limit []l9| . 

The condition that Aq ^ 1 should be more properly formulated as Aq <C 1/4, since A = 1/4 
is the middle point between the limiting saddle points at A = and A = 1/2. This criterion 
gives the following range of validity for the DLA solution ( |3.13| ) in momentum space 

16q^t < ln(A;i/A2). (3.14) 



(B) When r/agT is small, the saddle point is very close to A = 1/2 : 
1 . . r 



^0 



2 ' l5asT 



(3.15) 



which gives: 

I ~ f,^^sr \r ^2^, (3.16) 

where (5 = (-2^(1/2)) = 28C(3) and lo = 2V'(1) - 2^/^(1/2) = 41n2. 

When r = \n{k\/ h?), eq. (|3.16| ) corresponds to the usual solution to the momentum BFKL 
equation. On the other hand, when r = ln(rj^/^^), it also gives the solution of the linearized BK 
equation in coordinate space: 




,2 g(^a,r \-\n^{r\/e 



Mr{rA_)^\ -^ , „ exp<^ \n )■ (3.17) 
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The saddle point in eq. ( p. 15 ) remains close to 1/2 for all r such that 5 <^ 1/4. This is 
realized when (in momentum space, for definiteness) 



(3.18) 



which, as we shall soon discover, leaves a substantial window for genuine BFKL behaviour (in 
the sense of eq. ( p. 16 )) above the saturation scale. 

Since there is no overlapping between the conditions ( 3.14| ) and ( 3.18| ), we conclude that in 
the intermediate range where 8 ^ rjoigT ^ 16, the saddle point is relatively close to A = 1/4. 
Indeed, by using the approximation (3.11), one can estimate that Aq = 1/4 for r/agT ~ 14. 

Eq. ( p. 16 ) is the solution that is usually considered in applications of the BFKL equation at 
very high energy [^4|. It shows an exponential growth with r and "infrared diffusion" towards 
small k±_ momenta, which would become problematic in the formal high energy limit r — > oo 
at fixed r (which is, of course, outside the kinematical range considered here; recall eq. ( |3.3D .) 
But at sufficiently high energy, Q1{t) > Q^, and the BFKL equation gets supplanted by the BK 
equation (or the more general equations in Refs. |15, 17, ^), in which non-linear effects provide 
a natural solution to the difficulties of eq. (p.l6|) alluded to above 

(C) For r/oisT negative and large, the saddle point is close to A = 1: 



Ao ~ 1 - 5i 



5^ = 



(3.19) 



and the integral is estimated as 
/ ~ exp 



{r + 2^a,r(-r)}y^ 



<5? 



(3.20) 



This saddle point is interesting only in coordinate space, where r = ln(rj^/^^) < 0, and gives 
the dominant behaviour Mr{r±) at short distances: 



^frirx) = rj^A exp < 2JasTln' 



'A2 



ln(l/r2 A2) 



(3.21) 



where we have taken l"^ = 1/ h?. As expected, this is the same as the asymptotic solution of the 
DLA equation in coordinate space^ The saddle point ( 3.19 ) is close to 1 provided 6i <^ 1/4, 
which gives 



2 \2\ 



16asT « ln(l/rlA 



(3.22) 



in complete agreement with the previous condition (3.14) in momentum space. 
It is convenient to summarize the previous results in terms of the ratio 



R 



i-r) 



J-ln^ 

OoT A2 ' 



(3.23) 



^The additional factor r^A^ in front of the exponential in eq. (3.21) as compared to the DLA solution ( |j.l3[ ) 
in momentum space comes up since, in coordinate space, the DLA equation applies to Afr{r±) /r^. 
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(i? is positive in the perturbative regime of interest) which will play a special role in the next 
subsection. Our analysis shows that the saddle point solution to the BFKL equation has the 
genuine BFKL behaviour ( |3.16| ) for R <^ 8, the DLA behaviour (3.21) for R ^ 16, and some 



intermediate behaviour in the window 8 ^ ^ 16. These conclusions, based on a pure BFKL 
analysis, are still to be combined with the constraint ( |3.3| ) showing that we are indeed in the 
linear regime. This requires to determine the coefficient c in eq. (^.j), which we shall do in the 
next subsection. 

3.2 Saturation scale and scaling from the BFKL equation 

As already explained, it is possible to determine the saturation momentum from the solution 
■^T{r±) to the BFKL equation (|3.lD by using the saturation criterion: 



K{r± = 1/Qs(r)) = 1. (3.24) 

Strictly speaking, the number in the r.h.s. of this equation is not exactly 1 (cf. the discussion at 
the beginning of Sect. 3), but this is irrelevant to the accuracy of the present calculation. Indeed, 
replacing 1 by k < 1 would modify the following results via subleading terms, of relative order 
{l/asT)ln K,. Besides, such terms would affect the overall normalization of the saturation scale 
( |3.29| ), that is, the value of the reference scale A^, which has not been fully specified anyway. 
In what follows, we shall rely on the saddle point approximation in eqs. ( |3.8| )-( 3ri0| ) to obtain 



an estimate for Qs{t) 0, and then study the scaling properties of the BFKL solution above the 
saturation scale. For these purposes, it is crucial to notice from eq. ( |3.10| ), which is rewritten as 

^ U{\) + V;(l - A)) = ^ ^ -i? , (3.25) 

OA V /An a.T 



that the saddle point Aq is actually a function of only one variable R (cf. eq. ( 3.23 )) : 

Ao(r,r) = Ao(i?). (3.26) 

If one estimates the Mellin integral just by the saddle point, one then obtains 

Mr{ri_) ~ q^stF{\„{r),r) ^ ^3^27) 

where, as compared to eqs. (|3.8|) -( ^T0|) , we have changed the definition of the function F(A, r, r) 
by pulling out a factor a^r. This is convenient since, when evaluated at the saddle point, the 
new function F(Xq(R), R) is only a function of R. Also, we have neglected the (slowly varying) 
factor coming from the Gaussian integral over the fluctuations. This is correct up to corrections 
of relative order ln{asT)/asT. The effect of this factor will be illustrated in eq. ( ^.431) below. 
The saturation criterion ( 3.24D yields then the following condition on F : 



F{Xo{Rs),Rs) = for Rs = -J-ln^^ (3.28) 

which is an equation for Rs, and ultimately for Qs{t). 

This equation has two immediate and important consequences, which are among the main 
results in this paper: (i) The saturation momentum is increasing exponentially with r, with the 
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slope of the exponential uniquely fixed by the saddle point solution to the BFKL equation, (ii) 
The "anomalous dimension" which characterizes the approach towards saturation (cf. eq. (|1.3D ), 
is constant, i.e., independent of r, which then implies geometric scaling. 

i) Indeed, the solution Rg to eq. ( |3.28D is a pure number, Rs = c, and not a function of r 
(as it could have been if the function F{X, r, r) in eq. ( p. 27] ) was a function of two independent 
variables, r and r, and not just of their ratio R). Together with the second equation ( |3.28| ), this 
implie^: 

Q2(r) = A2e'=°=^ (3.29) 

which is consistent with a previous result ( |2.13| ) based on scaling properties of the BK equation, 
but which here arises from a (BFKL) solution which has no geometric scaling (e.g., eqs. ( |3.16| ) 



or (3.21) are not scaling functions for generic values of r± and r). Moreover, the value of the 



slope parameter c is known here, since uniquely fixed by the solution to eqs. ( 3.28| ) and ( p. 25 ). 



In fact, in order to compute Rg = c, there is no need to solve the saddle point equation ( |3.25| ) 
for arbitrary values of R (although this could be done via numerical techniques). Rather, we 
notice that by combining eqs. (1331 ) and (|^) one can deduce an equation for the particular 
value Xs = Xo{Rs) that the saddle point takes at saturation: 

- A.^(V(A) + V(l - A)) = XsRs = 2V'(1) - V(As) - ^'(l - A.). (3.30) 

We have solved this numerically and obtained: 

Xs = 0.627549..., R^ = c = 4.88339.... (3.31) 

Note that this Xg is not too far away from 1/2, which is consistent with the fact that the above 
value of Rs is in the range where we expect a "genuine" BFKL behaviour, i.e., Rs < 8. We shall 
return to this point after the discussion of geometric scaling. 

ii) Let us evaluate eq. ( |3.27| ) for R slightly above Rs, that is, for distances r± which, while still 
being much shorter than the saturation length 1/Qs{t), are nevertheless close to it in logarithmic 
units. More precisely, we shall require that (with Q'^ = l/r'j_) : 

0<R-Rs<^Rs, or l< In < In -^1^ . (3.32) 

The condition R > Rs (i.e., ^ (''")) ensures that we stay in the linear regime; in logarithmic 
units, this is effectively implemented as > eQl{T), with e = 2.72.... The condition R — Rs <^ 
Rs allows us to study the approach of Mr{r±) towards saturation in a limited expansion in 
powers of R — Rg. To linear order in this expansion, one has 



F{Xo{R),R) ~ F{Xo{Rs),Rs)+ ^F{Xo{R),R) 



{R-Rs) + --- 

R=Rs 

Xs{R-Rs) + ---, (3.33) 



^If we had included the Gaussian fluctuati ons a round the saddle point, this pure exponential would get mul- 



tiplied by a slowly varying prefactor; see eq. (3.43) for an example 
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where we have used the saturation condition ( 3.28| ), together with the fact that Ao(-R) is a 
solution of the saddle point equation, cf. eq. ( 3.10D , so that: 



dR 



F{Xo{R),R) 



R—Rs 



-HRs) + 



dF 
9A 



A=Ao 



R—Rs 



(3.34) 



Therefore, above the saturation scale, the dipole-hadron scattering amplitude is approximated 
as follows [recall that 6lst{R - Rs) = \ri{\ / r\Ql{T))] 



(3.35) 



where we have reintroduced the numerical factor k < 1 (cf. the discussion after eq. ( |3.24| )). 

Eq. ( |3.35| ) shows geometric scaling with the anomalous dimension^ 7 = 1 — A^ ~ 0.37 
determined by the value of the saddle point As = Ao(-Rs) at saturation. We summarize here the 
two essential ingredients in the arguments leading to this scaling: (a) the saturation condition 
which requires that Mr{r±) — > 1 when r± — > 1/Qs(t); (b) the fact that, at saturation, the saddle 



point Xo{R) becomes independent of r (i.e., a pure number; cf. eq. ( ^.31|) ). In turn, point (6) is 
the consequence of the "scaling" proper ty (|3.26| ) of the saddle point, and ultimately reflects the 
scale invariance of the BFKL equation (|3.l| ). 

The previous arguments also shed light on potential sources of scaling violation when going 
beyond the present approximations. For instance, we expect r-dependent corrections to the 
anomalous dimension due to the Gaussian fluctuations, the higher-order terms in the saddle 
point expansion, the initial condition Xo(A) in the Mellin function ( |3.6| ), etc. These corrections 
are suppressed by l/agT. Also, there are corrections to the functional form J\fT{r±), coming 
from higher order terms in the expansion ( p.33| ) around the saturation scale. These corrections 



are controlled by the ratio {R — Rs)/Rs = In^^/ln— ^ 
the second order in this expansion, eq. ( |3.35| ) gets replaced by 



For instance, after also including 



2asT 



In- 



1 



(3.36) 



where A'^ = {dXo{R)/dR)\R^. Clearly, the exponential term in this expression violates scaling. 
The scaling behaviour ( 3.35| ) holds in a window in specified by eq. ( ^.32| ), that is: 



1 < In 



Qi{r) 



CUsT <^ casT. 



(3.37) 



^We refer to 1 — As as an "anomalous dimension" since, naively, one could expect AfT{r±) to vanish like 
as r± 0: this would follow from the fact that A/'t(0) = 0, together with analyticity near r± 0. In reality, 
this analyticity is broken by logarithmic ultraviolet divergences in the formal expansion in powers of r^. These 
divergences correspond to the large logarithms ln((5^/A?) whose resummation by the DGLAP (or DLA) equation 
leads to the non-analytic behaviour manifest in eq. ( 3.21 ). Besides, whatever are the analytic properties of the 
function A/'-r(rx) near r± ^ (i.e., at very high the re is n o rea son why these properties should persist down 
to the = Qb(t). Compare, in this respect, eqs. (3.35) and (3.21). 
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This requires that both r = ln(Q^/A^) and cagT are large numbers, with r larger than cagT (in 
agreement with the original assumptions (3^) and ( |3.3D ), but not much larger: r — cagT <^ cagT. 
In practice, this requires ln((5^/A^) to be deeply inside the strip (note that c = 4.88... ~ 5) 



5a,r < In(QVA^) < WosT , 



(3.38) 



which is indeed a large window when Sa^r » 1. Together with eq. ( 3.18| ), this suggests that 
the window for extended scaling is almost entirely located in the kinematical range controlled 
by the BFKL saddle point ( jslsl) . 

To verify that, let us compute directly the predictions of this saddle point and the associated 
solution ( ^.17 ) for the saturation scale and geometric scaling. The "BFKL saturation scale" is 
defined by imposing the condition ( |3.24| ) directly on eq. ( p.l7| ): 




1. 



(3.39) 



As expected, this amounts to an equation for Rg, the value of the variable R at saturation 
(compare to eq. ( p. 28 )). In this case, this is a second-order equation: 

0, (3.40) 



R; + PRs - 2f3uj 



P + ^I3{P + 8uj) 



4.8473. 



(3.41) 



with the positive solution: 
-Ro = — 

BFKL 2 

This is the value of the slope parameter c predicted by the standard BFKL solution, and is 
numerically very close to that in eq. ( p. 31 ). By using ( |3.41| ) and ( p. 15 ), one can compute the 
BFKL saddle point at saturation As (or the anomalous dimension 7 = 1 — As): 



A< 



BFKL 



1 -R. 

2 /5 



0.644..., 



(3.42) 



which is indeed very close to the true saddle point ( |3.3lD . This shows that for all practical 
purposes one can use the explicit BFKL solution ( |3.17| ) for any = l/r'j_ in the window for 
extended scaling. In fact, without any approximation, eq. ( p. 17 ) can be cast in the form of the 
"second-order expansion" in eq. (|3^) , with As given by eq. (|3.42|) , and A^ = 1//3 ~ 1/33.67. 
This latter is a rather small number, showing that the violations of scaling due to the exponential 
term are only tiny. That is, in its whole domain of applicability, the "genuine" BFKL solution 
( 3.17| ) is almost an exact scaling solution. 

For completeness, and also for comparison with previous analytic studies in the literature 
which have used the DLA approximation |8|, ^, |l^] , let us finally evaluate the scaling predictions 
of the DLA saddle point ( 3.19 ). The saturation condition ( 3.24 ) applied to eq. ( |3.21 ) yields: 



A^ 



(3.43) 



DLA 327rQ!sT 

where the slowly varying factor multiplying the exponential is the effect of the Gaussian fluctu- 
ations around the saddle point. Eq. ( |3.43| ) implies: 



R. 



ln(327rasT) 



DLA 



4, 



(3.44) 
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which is slowly dependent upon r, because of the contribution of the fluctuations, but to leading 
order takes the constant value cdla = 4. This is in agreement with previous studies [^, ^ 
If inserted in eq. (3.19), it gives the DLA anomalous dimension 7 = 1 



A., 



7 



DLA 




1^1 ln(327rasr) 



(3.45) 



This is slightly larger than for the true saddle point, or the BFKL saddle point ( 3.42| ). The DLA 
window for extended geometric scaling is given by eq. ( 3.37 ) with c = 4. 

Thus, at a first sight, the DLA predictions look rather similar to those of the BFKL approx- 
imation, both for the value of the saturation scale (which fixes also the scaling window), and 
for the anomalous dimension. However, one should keep in mind that these predictions are not 
consistent with the validity region of DLA: the saddle point ( 3.45| ) not only is not close to 1, 
but it even takes the typical BFKL value. Besides, the corresponding scaling window is truly 
within the kinematical range for BFKL. 

We finally mention that numerical studies of the BK equation |^, |ll| (with a fixed 
coupling) have found that the saturation scale is indeed increasing exponentially with r, with a 
slope parameter c ~ 4.1 which is intermediate between the DLA and BFKL predictions obtained 
in this section. 



3.3 The running coupling case: ^^((^^(t)) 

Except for the modified r-dependence of the saturation scale, which changes according to the 
general expectation in eq. ( p.20|) , all the previous discussion of the solutions to the BFKL 
equation and their consequences for extended scaling goes almost unchanged to the case of a 
running coupling in which the scale for running is set by Qs ("?")• In this subsection, we shall only 
indicate the few steps which involve non-trivial differences. 

The extra r dependence of the coupling as{Q'^{T)) commutes with the Mellin transform 
( |3.4| ), which is defined only in terms of the transverse coordinates. Thus, the previous equation 
for Xt(A) is simply replaced by 

l-^-(^) = r77527T7A^^2^W-^(A)-V'(l-A)}xr(A). (3.46) 
dr ln(Qf(r)/A^CD) 

So far, the saturation scale is not known, but it can be absorbed into a redefinition of the "time" 
variable r : 

ln(Q?(r)/A|cD) ^ - ^, or f ^ £ dr' , (3.47) 

where we have written QI{t) = Aqq^b-'^^^^ Then, the equation is easily solved: 

X.(A) = e^of {2^(i)-V'(A)-V(i-A)} ^ (3^48) 

which leads us to the following Mellin representation for the solution to the BFKL equation 
with running coupling: 

Mr(r±)= e^(^'''^) with F(\,r,f) = rX + bof{2Ml)-'il^(X)-m-X)} (3.49) 

Jc 2m 
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where r = ln(r^/£^) = In(r^A^) as before. As anticipated, this has the same structure as in the 
fixed couphng case (cf. eq ( |3.7D ). Thus, all the results in Sects. 3.1 and 3.2 can be immediately 
translated to the case of a running coupling by simply replacing 



and 



(3.50) 



in the corresponding formulae. In particular, by the same arguments as before, the saddle point 
Ao is a function of the ratio R alone, with R = {—r)/bQf. The value = Xo{Rs) of the saddle 
point at saturation is again determined by eq. ( 3.3C| ), so that A^ and c = Rs take the same 
values as before, cf. eq. ( |3.31| ). Thus, all the previous results on extended scaling — the value 
of the anomalous dimension 7 = 1 — A^ in eq. (|]35D, and the momentum range (|3.32| ) in which 
the scaling holds — remain unchanged, except for the expression of the saturation scale which 
enters these results. 

To determine this scale, we use c = Rs = (— rs)/6of with (— r^) = lnQ^{T)/ Aqq^^ = /(r). 



together with eq (|3.47D , to derive an equation for /(r) 



1 



fir) 



bo f dr'-^. 
Jo f[r') 



This has the solution /(r) = A/26oc(r + tq) (and therefore f 
saturation scale is determined as 



Qs(^) = Aocoe 



(3.51) 

\/2{t + To)/bQc). Therefore, the 



(3.52) 



where tq is arbitrary and c is given by eq. (3.31). This has the same functional form as eq. ( |2.20| ) 
that has been obtained from the BK equation with a scaling Ansatz. In particular, this confirms 
the previous arguments in Sect. 2 that the coefficient c which enters the exponent of the 
saturation scale should be the same for fixed coupling and running coupling. 

Needless to say, the BFKL and DLA predictions for the anomalous dimension, eqs. ( |3.42| ) 
and respectively ( 3.45|) , remain unchanged. 



4 Summary and discussion 

We have shown that the geometric scaling predicted at low momenta ^ by the Colour 
Glass Condensate and phenomenological saturation models is preserved by the BFKL evolution 
equation up to relatively large momenta, within the range 1 ^ ln(Q^/Qg) <C ln{Q'^ / Aqq-^) . 
By matching the solution to the BFKL equation onto the saturation condition at '-^ Q^, we 
have determined the dependence of the saturation scale Qs upon the rapidity r = ln(l/x), and 
the anomalous dimension of the distribution function near saturation. 

The matching has been performed for the dipole scattering amplitude, which enters linearly 
the structure function F2{x,Q'^), and can be also related to the gluon distribution. We have found 
the same kinematical window for extended scaling for both fixed and running couplings (with the 
scale for running set by Qs{t)), although the r-dependence of the saturation momentum turns 
out to be different in the two cases. We have also found that, formally, the double logarithmic 
approximation to the DGLAP equation leads to qualitatively, and even quantitatively, similar 
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predictions for the extended scahng, but these results are inconsistent with the validity range of 
this approximation. 

We have shown that the functional form of the r-dependence of the saturation scale Qs 
which follows from the BFKL equation is consistent with the general scaling properties of the 
BK equation. 

It would be extremely interesting to compare these results with the F2 data in deep inelastic 
scattering, for which geometric scaling has been originally identified Q, and to explore possible 
implications for particle production in heavy-ion experiments, where a phenomenological scaling 
law has been recently reported |35]. In particular, we expect these results to have consequences 
for the analysis of the multiplicity distributions of produced particles in heavy ion collisions 
at RHIC Our analysis shows that while F2 has geometric scaling, the gluon distribution 
function, due to an extra factor of l/ag in eqs. ( p.7[ )~([2^), has logarithmic scaling violations. 
This is consistent with scaling violations which are seen in the RHIC data 

While the qualitative phenomenon of the existence of extended geometric scaling up to mo- 
menta of order 100 GeV^ is certainly consistent with the analysis of deep inelastic scattering 
by Stasto, Golec-Biernat and Kwiecihski Q, it is at the same time clear that the saturation 
scale emerging from our analysis (for either fixed or running coupling) is too rapidly increasing 
with T = ln(l/x) to give a good description of the data for F2. 

Consider fixed coupling first. Although eq. ( p. 29 ) shows a power law increase, Ql{x) = 
Q'qX~^, in agreement with the parametrization used in Ref. H], the actual value of the parameter 
A predicted by the BFKL equation, namely A ~ (4 — h)asNc/'K, is sensibly larger (for Nc = 2> 
and realistic values of as) than the value A = 0.3 — 0.4 extracted from the fit to the data [^. 

As for the corresponding prediction in the case of a running coupling, eq. ( p. 52 ), this is less 
rapidly increasing at very large r : Q1{t) = ^^q-q exp{-\/C(T + tq)}. But since the number 
C = 2boc is relatively large (of order 10; cf. eqs. ( p.l8| ) and ( |3.31| )), this too fails to reproduce 
the data, unless the parameter tq is mysteriously high. 

The reasons for such a failure can be several. First, there are the many approximations 
that we have performed in order to obtain analytic solutions to the BFKL equation. Note that 
we have preserved just leading order terms in the exponents, and that the subleading terms 
are truly small only if the inequalities ( |3.2| ) are strictly satisfied. These are strong inequalities 
on logarithms, which may not be well fulfilled for realistic values of x and Q^. Besides, even 
small subleading terms may give a substantial effect once exponentiated. This is illustrated by 
eqs. (|3.43D -( ^.45|) which, in addition to the lowest order terms, include also the corrections due 
to the Gaussian fluctuations of the saddle point. As obvious on eq. ( |3.43| ), these "subleading" 
terms may drastically change the actual value of the saturation momentum. Still, the fact that 
the slope parameter c ~ 4 — 5 that we have found is rather close to that obtained via numerical 
studies of the BK equation [10,11] makes us confident about our control of the approximation 
scheme for the exponent. 

The last argument also suggests that, independent of further approximations, there is a true 
discrepancy between the exponent A ~ 4a<j predicted by the BK equation and its phenomeno- 
logical value A = 0.3 — 0.4 [^, 23|. Recall that the BK equation is strictly valid only in the 



large Nc limit, and to lowest order in as (within a leading- log approximation scheme). For finite 
Nc, there is no simple closed non-linear equation, just an infinite hierarchy of coupled equations 
1 15 1, which is equivalent to a functional Fokker-Planck equation ^. This functional equation 
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can be studied via numerical techniques (in particular, on the basis of the associated Langevin 
equation |17|, or on the path integral representation |36|), and it would be very interesting to 
estimate the finite Nc corrections. Moreover, higher-order corrections in Og, which so far have 
not been included in the non- linear evolution equations (see however [^]), but which are now 
available for the BFKL kernel [^], may be responsible for an effective decrease in the slope 
parameter A with respect to its lowest order BFKL value. 

Note finally that the uncertainty on the r-dependence of the saturation scale should not 
affect our prediction for the kinematical window in which one expects extended scaling. Indeed, 
as explained in relation with eq. (|1.4| ), this prediction relies just on a limited expansion around 

— Q% together with the scale- invariance of the linear evolution equation at hand (in our case, 
BFKL). This argument predicts an upper limit Q^ax ~ Qs/^qcd which geometric scaling 

should be expected. To estimate this upper limit, we shall use phenomenologically reasonable 



values for Qs [23, ^, and not the theoretical predictions of our analysis (which cannot give 
the absolute value of Qs anyway, just its r-dependence). This gives Qmax ~ 100 GeV^ for 
Qs ~ 1 GeV, and Q^„^ ~ 400 GeV^ for ~ 2 GeV, values which are both of the right order 
of magnitude to agree with the phenomenological analysis in Ref. [0]. 
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Note added 

Very recently, after this work was essentially completed, a paper has been released by Kwiecihski 
and Stasto which addresses the issue of geometric scaling at high Q^ in the framework of 
the DGLAP evolution equation (or its double-log approximation). The conclusions in this paper 
are however different from ours: The authors of Ref. have not recognized the existence of a 
window for geometric scaling above Qs, but rather concluded that scaling violations should be 
expected for any Q^ > Q'^{x), even at very small values of x. 
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